Non-fermi liquid behavior in itinerant antiferromagnets 
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We consider a two dimensional itinerant antiferromagnet near a quantum critical point. We show 
that, contrary to conventional wisdom, fermionic excitations in the ordered state are not the usual 
Fermi liquid quasiparticles. Instead, down to very low frequencies, the fermionic self energy varies 
as to 2 / 3 . This non- Fermi liquid behavior originates in the coupling of fermions to the longitudinal 
spin susceptibility x\\(q,ty m which the order-induced "gap" in the spectrum at q = dissolves 
into the Landau damping term at v F q > Q. The transverse spin fluctuations obey z = 1 scaling 
characteristic of spin waves, but remain overdamped in a finite range near the critical point. 



Fermi liquid (FL) theory is a powerful tool to study 
properties of interacting electrons. It states that, upon 
switching on the electron-electron interaction, elemen- 
tary excitations near the Fermi surface remain fermionic, 
with well-defined energy and momentum. As a result, 
Fermi liquids have a constant uniform spin susceptibility 
and a linear in temperature, T, electronic specific heat. 

In the last two decades, many new compounds were 
found which exhibit thermodynamic and transport prop- 
erties dramatically different from those of a Fermi liq- 
uid . Understanding the origin of this non- Fermi liquid 
(NFL) behavior is one of the most important challenges 
in condensed matter physics. 

Few avenues leading to the NFL behavior in dimen- 
sion D > 1 have been proposed in the theory of in- 
teracting electron systems. One of the most widely 
studied scenarios ties the destruction of the Fermi liq- 
uid to a quantum critical point (QCP) @. At a QCP 
the fluctuations of the order parameter(OP) are gap- 
less, and the effective electron-electron interaction me- 
diated by this bosonic mode is long-ranged. Scattering 
of electrons by gapless bosons is singular and destroys a 
Fermi liquid when the local static bosonic susceptibility 
Xi= J d D - 1 q x{q,u = 0) diverges 0. 

It is generally expected that, away from the QCP, the 
Fermi liquid behavior is preserved. On the disordered 
side the FL is protected by the gap in the OP fluctuations 
due to a finite correlation length. On the ordered side, 
where the OP acquires a mean value, the amplitude (lon- 
gitudinal) bosonic excitations are gapped, while the gap- 
less Goldstone excitations (phase, or transverse modes) 
are harmless for fermions due to the Adler's principle |4j 
that states that the fermion-boson vertex vanishes at the 
momentum transfer equal to the ordering wave vector Q. 

The main conclusion of this Letter is that, contrary 
to this general belief, a novel non-Fermi liquid electron 
behavior emerges on the ordered side of a QCP when Q ^ 
0. This NFL behavior originates from the interaction 
between electrons and longitudinal bosonic excitations. 
The key observation is that this mode becomes gapless 



and overdamped in the range v F q > ^ due to Landau 
damping from the electrons, and gives rise to a strong 
electron-electron interaction even at the lowest energies. 
The NFL behavior persists down to a frequency uj m i n 
that vanishes if the fermionic bandwidth is infinite. 

Below we consider a SDW transition between a para- 
magnet and an itinerant antiferromagnet (AFM) with 
Q = (-7T, 7r) and the dynamical exponent z = 2 in D = 2. 
Such a transition is both one of the most studied and the 
most relevant experimentally, as it occurs in many heavy 
fermion materials and is believed to be responsible for 
their unusual properties. 

We first present the results and discuss the physics 
and then report the details of the calculations. We mea- 
sure the AFM order by the gap in the fermionic spec- 
trum, to. Opening of this gap renormalizes the spin sus- 
ceptibilities from their form on the paramagnetic side, 
Xl^q.fi) = xf(q,n) = (q - Q) 2 ~ 1- Here- 

after q denotes the deviation from (n, 7r), except when 
noted. Conventional wisdom holds that, in the AFM 
state, at energies smaller than 2to, the transverse spin 
excitations are the Goldstone spin waves, while the lon- 
gitudinal excitations are gapped at all bosonic q and Q. 
We find, however, that this behavior holds only when 
v F q is smaller than the frequency, 2to > Q > v F q. In 
this range, xj ~ q 2 — a^V? jm, has a spin- wave form (a 
is a number O(l)), while X\\ X ~ I 2 + a^^Jm 2 — f2 2 /4 has 
a gap, 2m, for the propagating spin excitations. 

In the opposite limit, 2m > v F q > f2, we find 

Xjj" 1 ~ q 2 -2i-fm\n\/v F q, Xj 1 ~ <? 2 - ij\tt\v F q/m. (I) 

The transverse spin excitations become overdamped, 
though they still show z = I scaling. More importantly, 
the constant 2m gap in the longitudinal spin excitations 
dissolves into the Landau damping term, and the longi- 
tudinal excitations become gapless with the dynamical 
exponent z — 3. This implies that at a given v F q <C to, 
there are two regions of (f2 < v F q and £1 > 2m — v F q), 
where x'\\{Qi ^) is nonzero (see Fig- EJ- 
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FIG. 1: Doubling of the unit cell in an AFM transition. Verti- 
cal dotted lines are the boundaries of the magnetic Brillouin 
Zone. The Fermi surface at the "hot spots" at the para- 
magnetic side of the transition is the dashed red line, the 
dash-dotted blue line is its translation by Q. Solid black lines 
denote the Fermi surface in the AFM state. 

The gapless overdamped form of the longitudinal sus- 
ceptibility strongly affects fermionic self-energy as the 
interaction vertex for x\\ ls n ot reduced by Adler's sym- 
metry principle. Critical theories with z = 3 have been 
studied extensively [|J. In two dimensions, the electron 
self energy due to interaction with such a mode varies as 
(iuj) 2 / 3 , i.e. the Fermi-liquid behavior is violated despite 
the presence of the AFM order. This self-energy comes 
from vpq oc SI 1 / 3 ^ SI where Eq. is valid. Similar 
power law for the self energy in the ordered state was 
argued to be present in electronic nematic phases ||. 

Above results are for an infinite fermionic bandwidth, 
W. When W is finite, the longitudinal spin fluctuations 
have a gap of the order of m 2 /W even for SI < vpq. 
This scale sets the lower cutoff for lu 2 ^ 3 behavior; the 
Fermi liquid behavior is restored at lower energies. It is 
essential that in itinerant AFM, m <C W (see below), 
and the cutoff scale is parametrically smaller than m. 

The evolution of the spin response with tt/vpq is re- 
lated to the fact that the mass term in x\\ (<Z; ^) is propor- 
tional to the particle-hole polarization bubble H(q, SI) = 
II(g, SI) — II(q, 0), see below. In a nested antiferromagnet, 
the AFM order opens a gap over the entire Fermi sur- 
face. In this case, IT(q, SI) is a constant independent on 
the ratio Sl/ii^g. However, in an itinerant AFM without 
nesting, only parts of the Fermi surface close to the "hot 
spots" become gapped due to long-range spin ordering. 
Elsewhere on the Fermi surface a continuous particle-hole 
spectrum still exists. As a result, spin excitations can de- 
cay into particle-hole pairs. On the paramagnetic side, 
this process involves large momenta^ near Q, hence the 
damping is linear in frequency, ImIT(q, SI) = 7|S7|. In 
the ordered state, the lattice period is doubled, and the 
fermionic momenta k and k + Q become equivalent. As 
a result, a bosonic mode at the momentum Qi ~ Q scat- 
ters a fermion from a point kp into a Fermi surface point 
k^+q with small q = Qi Q, see Fig. B The polarization 
bubble then becomes n(q, SI) oc 7771SI/ a/S! 2 — (vpq) 2 . It 
gives rise to a mass term for S7 3> vpq, but reduces to 
ijmfl/vpq in the limit SI <C vpq. 



We now provide the details of the calculation. To in- 
vestigate the properties of an itinerant antiferromagnet 
near a QCP we employ the spin fermion model, which 
has been widely used to study the properties of param- 
agnetic metals close to a magnetic instability 0- Here 
we use it on the ordered side of the transition. The model 
Hamiltonian consists of three parts: the electronic band, 
H^p = Eka £ W c ka c M' the collective fermionic spin 

excitations, — X) q Xo 1 (q)S q S_ q , and the interac- 

tion of electrons with the spin fluctuations, 

H sf = 9 Clcc^SqCk+q^ + h.C. (2) 

k,q,a,/3 

Here e(k) is the bare quasiparticle dispersion, xo(q) is 
the bare static spin susceptibility and g is a coupling 
constant that we take to be independent of momentum. 
The rationale for this model is that the antiferromagnetic 
interactions emerge from the energies up to the fermion 
bandwidth. Therefore to analyze the behavior of low- 
energy quasiparticles it is sufficient to separate the energy 
scales and treat the collective spin degrees of freedom 
as a separate bosonic field. This procedure is internally 
consistent if the coupling g <^W. 

Properties of the spin subsystem are determined by the 
bare susceptibility, which we take to have the Ornstein- 
Zernike form. xo(q) = Xo/(£q + (q _ Q) 2 )- as in earlier 
work In the paramagnetic phase, the magnetic corre- 
lation length, £q, is real and controls the separation from 
the AFM instability. The full dynamical susceptibility is 

X(q. w) = Xo/tto 2 + (q - Q) 2 + H(q, w)), where n(q, w) 
is the full polarization bubble that has to be evaluated 
within the low-energy theory. 

On the antiferromagnetic side of the transition £ be- 
comes imaginary < 0), so that the susceptibility 
diverges near the AFM wave vector, Q, leading to a 
staggered magnetic moment, (Sq). The divergence is 
compensated by n(q, 0) which acquires a finite value in 
the ordered phase. Furthermore, for a finite (Sq), the 
longitudinal (zz) and the transverse polarization bub- 
bles differ, leading to the anisotropic susceptibility. The 
staggered moment, (Sq), is related to £q~ 2 < 0, via 
the Goldstone requirement that the fully renormalized 
xl 1 ^ 0,0) = 0. 

Following the standard procedure, we introduce m = 
(/(Sq), include the "condensed" part of the spin fluctu- 
ations, ^ k a ocmc^Q Q Ck.ai into the fermion part of the 
hamiltonian, and diagonalizc it by Bogoliubov transfor- 
mation to find H f = J2n,k,a ^n( k ) a i,k,Q a ™,k,Q- The dis- 
persion of the fermions in the two new bands is given by 
£±i(k) = e k ± ^S 2 + m 2 with e k = (e(k) + e(k + Q))/2 
and <5k = (e(k) — e(k + Q))/2. The new and the original 
fermion operators are related by a unitary transformation 
a li.k. Q = "kcL - <WkC k+QiQ , and al i k a = av * c l,a + 
«k4+Q,a. with ul(v 2 ) = [1 + (-)S k /^/S 2 +m 2 }/2. The 
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interaction of fermions with the uncondensed part of the spin excitations now takes the form 
I 



H sf =9 T nf,n 2 ( k , k + q)ai l!k!Q an 2 ,k+q,/30-a/3S q _Q + h.C. 

ni ,ri2 ,k,q,a./3 



r^„ 2 (k,k + q ) = 



MkUk+q - a/?fkWk+q 



1 fini,ri2 



(3) 
(4) 



The vertex for the Goldstone transverse mode, r" : "(k, k + q) cx |q| at small q, as required by Adler's principle. The 
spin polarization bubble U(q, Q) is given by 



where Gj(k, u) n ) is the Green's function of the fermion 
in band i at the Matsubara frequency, u> n = irT(2n + 
1). The tensor Hjj is diagonal, and its zz (xx and yy) 
component is the longitudinal (transverse) polarization. 

The change in the static polarization upon entering 
the AFM state, U(q -> 0,0) - n m=0 (g -> 0,0) extends 
over all energies up to the fermionic bandwidth. This is 
expected since this difference has to compensate negative 
£o~ 2 that also comes from high energy fermion physics. 
The relationship between m and Q 2 then depends on 
the regularization chosen for the high energies. We will 
simply measure the distance from the QCP on the AFM 
side in terms of m, and use the Goldstone condition £(f 2 + 
II_l(<? — * 0,0) = 0. In contrast, Hij(q, f2) are independent 
of the cutoff procedure and are fully accounted for in the 
low-energy theory. It is essential that Tl\\(q — * 0,0) — 



U±(q -> 0,0} = 0(1/W), and hence to this order, £ 2 + 

ny(9,n) = n u (g,n). 

Evaluation of ILi and II j_ is straightforward. We con- 
sider a two-dimensional Fermi surface with tetragonal 
symmetry, so that there are four pairs of hot spots which 
contribute additively to H(q, f2). Without loss of gen- 
erality we choose the Fermi velocities at the hot spots 
kp and kp + Q along x and y directions, respectively. 
It is convenient to introduce k± — (k x ± k y )/2, so that 
£k = v F k + and <5k = v F k-. These directions are inequiv- 
alent as the Fermi surface exists for arbitrary k_, but 
only for |k+| > m/v F> see Fig.l. Evaluating the inte- 
grals in Eq.© we find for a given pair of "hot spots" 

(n > o) 



IIx(q,Q) 



n.,(q,n) = 



r yflsgn(v F q + — Q) 



7ilsgn(57 — vpq + ) 



where 7 = g 2 xo/(^ 7TV F ), and Q = fl + iSn- The full 
II(g, Q) are the sums over the pairs of hot spots, i.e. the 
sum of contributions from the wave vectors (±q+ , ±<7- ) 
and (±<?_, ±<Z+)- 

The behavior of the susceptibilities is particularly sim- 
ple if we set q+ = in Eqs.®-(|ZI)- When Q,v F q_ > to, 
Uj_(q,n) = 8||(cj>, Q) ~ i'7|fi|, as at the QCP with 



= 2. For to > ft > v F q- 



n 



n 2 /m, 



i.e. 



Xj_ cx qi 



(7/771)0 . In the same regime ILi ss 27TO, 



Am 2 



\v F q+ - n) 2 - v 2 F q 2 _ 
4m 2 

(v F q + - Q) 2 - v 2 F q 2 _ 
I 



-1/2 



-1 1/2 



(0) 
(>) 



and therefore longitudinal excitations have a gap. In the 
opposite limit to v F q 3> O, we find II ^ oc i\Q\(v F q)/m 
and IX || cx im£l/v F q. Consequently, the two spin suscep- 
tibilities have the forms given in Eq. JQ, i.e., spin exci- 
tations are overdamped and gapless, with z = 1 {z = 3) 
for the transverse (longitudinal) channel. 

That undamped spin waves exist only for Q > v F q 
means that propagating undamped spin waves exist 
only above a threshold value of staggered magnetiza- 
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FIG. 2: Transverse spin susceptibility at £l,VFq <C m. The 
propagating spin wave is manifested as a delta-function peak 
at £l(q) in the top panel. 
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FIG. 3: Dissipative part of the longitudinal susceptibility. Be- 
low energies ~ 2m the onset of damping occurs at = VFq. 
Notice the break in the horizontal scale. Inset: low energy 
behavior of Imxn(q, fi). The maximum is at !) raal oc q 3 . No- 
tice that in the Landau-damped region Imx||(q, O) is much 
larger than that in the continuum. 



tion, m c ~ jVp ~ g (i.e., (Sq) = 0(1))- For m < m c 
the spin wave would be at Q < vfQ, but it is damped 
there. At m = m c the spin-wave pole first splits from 
the upper edge of the low-frequency continuum and a 
propagating spin- wave emerges, initially with the residue 
Z cx m — m m i n . This behavior is illustrated in Fig. [21 
For x\\i there is no distiction between different values of 
m - in all cases, for vpq <C m, there are two separate 
regions where dissipation is present: one is at f2 < vpq-, 
and the other at fl > 2m — Vpq—, as shown in Fig. |3| 

For a general q, when both q + and g_ are compara- 
ble, the details of the momentum and frequency depen- 
dence of the susceptibility are somewhat lengthly, but 
follow straightforwardly from Eqs.@-(71l. The Landau 
damping is now effective for fl < vp(\q + \ + \q~\) <C 2m. 

At small frequencies x^ 1 a? 2 - ijrnfl/(vp^J\q^_ — 
The threshold for the spin-wave creation along the direc- 



tion q + = aq- is m c = jvp(a + l) 3 ^ 2 (y/a + 1)/(1 + a 2 ), 
and the spin-wave pole first appears along the diagonals 
of the Brillouin Zone, at m > jVp. 

The regime Q <C vf|<7+ — Q-\ "C m is particularly im- 
portant for the electron self-energy, E(k, u>). To the low- 
est order in the interaction, 



x X i i (k-q ! a;-n)G„ 1 (q,fi). (8) 



We find that there are two non-trivial contributions 
to the self energy. The first comes from the large 
momentum transfers vpq > 2m and gives Ej 1 « 
(ig 2 \p2j 'vp)(y / 'Am 2 — ijauj — 2m) cx oj/m at small ener- 
gies. At QCP this term gave rise to the \fiuj dependence 
of the self energy at hot spots 3] . Finite staggered mag- 
netization, m, renders it harmless for the fermions. 

The second contribution appears for m =/= from 
the interaction between the electrons and the Landau 
damped longitudinal spin excitations. This contribu- 
tion is dominated by vpq cx J7 1 / 3 for which the cri- 
terium vpq > O is satisfied. In contrast to the trans- 
verse Goldstone mode the interaction vertex for X|| is 
not reduced by Adler's symmetry principle. As a re- 
sult, the corresponding contribution to the self-energy is 
E? w (7r/4v^HWMM 2/3 /( m 7) 1/3 , as in z = 3 crit- 
ical theories jfj. The fractional ui dependence with ex- 
ponent 2/3 implies that both 7?eEi and im£i are of the 
same order, giving rise to a non-Fermi liquid behavior. 
Furthermore, as the Landau damping is due to small mo- 
mentum transfer, the NFL form of the self-energy persists 
everywhere on the Fermi surface, not only at the former 
hot spots. This is not surprising since, in the presence 
of the long-range AFM order, the entire Fermi surface 
becomes "hot" as the points kp and ki? + Q become 
equivalent. Finally, in complete analogy to z — 3 gauge 
and ferromagnetic critical theories , vertex corrections 
and the momentum-dependent part of the self-energy are 
small for m < g <C W. Hence the result E cx (— zw) 2 / 3 
remains valid to all orders in the interaction. 

We believe that this scenario of Fermi liquid break- 
down on the ordered side of QCP is quite general. The 
requirement that the gap in the optical mode of the or- 
der parameter dissolves into the Landau damping due to 
small momentum scattering in the enlarged unit cell is 
satisfied not only for AFM, but also for CDW and other 
transitions with a finite Q. 

The most obvious experimentally observable conse- 
quence of this scenario is the existence of the two regions 
of £1 where Imx\\(<l, O) (see Fig. |3J). This can be 
probed by inelastic neutron scattering experiments. The 
fractional exponent in the self energy also leads to the 
non-Fermi liquid, quantum-critical behavior for the tem- 
perature dependence of the electronic specific heat and 
resistivity pj, which is accessible experimentally. 
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